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Abstract 

We quantize the recent hydrodynamic analysis of Stringari for the low-energy 
collective modes of a trapped Bose gas at T = 0. This is based on the time- 
dependent Gross-Pitaevskii equation, but omits the kinetic energy of the den- 
sity fluctuations. We diagonalize the hydrodynamic Hamiltonian in terms of 
the normal modes associated with the amplitude and phase of the inhomoge- 
neous Bose order parameter. These normal modes provide a convenient basis 
for calculating observable quantities. As applications, we calculate the deple- 
tion of the condensate at T = as well as the inelastic light-scattering cross 
section S(q, to) from low-energy condensate fluctuations. The latter involves 
a sum over all normal modes, with a weight proportional to the square of the 
q Fourier component of the density fluctuation associated with a given mode. 
Finally, we show how the Thomas-Fermi hydrodynamic description can be 
derived starting from the coupled Bogoliubov equations. 

PACS numbers: 03.75.Fi, 67.40.Db, 78.35.+C 
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I. INTRODUCTION 



The recent achievement of Bose-Einstein condensation (BEC) in trapped atomic gases 
at ultra-cold temperatures |IJ has stimulated interest in the collective oscillations of a non- 
uniform Bose condensate These are of direct experimental interest Q. At temperatures 
well below the BEC transition temperature, all atoms are in the condensate, and one may 
then base the analysis on the Gross-Pitaevskii (GP) equation of motion M for the time- 
dependent condensate wave function $(r, t) ||. Linearizing around the static equilibrium 
value $o( r )> this leads to the well-known Bogoliubov equations of motion for the excitations 
of the condensate (see, for example, Refs. j7|J§). Recently Stringari || has noted that the 
GP equation can be rewritten as equations for the condensate density and phase variables 
and that within a Thomas-Fermi approximation (TEA) which neglects the kinetic energy 
associated with density fluctuations, these have the familiar structure of the hydrodynamic 
equations of a superfluid at T = 0. This allows one to give a simple theory of the low- 
frequency excitations of a trapped Bose gas, without the necessity of solving the full set of 
coupled Bogoliubov equations 

In the present paper, we develop the approach of Ref. by quantizing this "hydro- 
dynamic" description and diagonalizing the associated Hamiltonian in terms of the normal 
modes of the condensate. This is a natural generalization of the well-known discussion used 



in a uniform Bose-condensed fluid [|T0]] . These modes form a natural basis for understanding 
the effect of low-energy condensate fluctuations on various physical quantities. We illustrate 
this in this paper by evaluating the local depletion of the condensate at T = as well as 
the inelastic light scattering cross section from the low-energy condensate fluctuations. The 
latter is found to be quite different from the case of a uniform weakly-interacting Bose gas 
(as discussed, for example, in Refs. |TI| , |12| ). 

The present "hydrodynamic" description is only an adequate description of the low- 
energy excitations of the condensate at T = 0. In particular, as discussed by Stringari 
0], the neglect of the kinetic energy in the TF approximation is only correct in the strong 
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interaction or large-density limit |L3| . In Section V, we discuss the precise relation between 
the present analysis and that based on the full coupled Bogoliubov equations for condensate 
excitations. In spite of its deficiencies, the hydrodynamic approximation has the advantage 
that it allows one to illustrate the qualitative effects of the low-frequency modes in a very 
explicit manner. 



II. QUANTIZED HYDRODYNAMIC THEORY 



In a dilute trapped Bose gas |]] near T = 0, most of the atoms are Bose-condensed and 
thus the whole system can be well described by the Gross-Pitaevskii Hamiltonian [|]|| 



H = J dr $*(r) 



$(r). 



(1) 



In ([[]), the static external potential V ex (r) is responsible for trapping of the atoms. We 
have assumed a short-range interaction g = A7rh 2 a/m between the atoms, with a being the 
s-wave scattering length. The Bose-condensate order parameter $(r) can also be written as 



i</>(r) 



$(r) = |$(r)|e 



in terms of amplitude and phase components. Substitution of (||) into ([[]) yields 



(2) 



H = dr 



h 2 \ 



T c 1 J V z ^p~ c + -mp c v 2 - p c ji + p c V cx + -g p 2 c 



where the local condensate density and superfluid velocity are defined by 



(3) 



v r 



2mi|$(r)| 2 
W0(r) 
m 



[$*(r)V$(r) - W(r)$(r)] 



(4a) 



(4b) 



Expressed in terms of these variables, the GP theory can be formulated in terms of 
hydrodynamic-type equations. We can also consider a time-dependent order parameter 
$(r, t) given by the time- dependent GP equation || 



h 2 v 2 



2m 



-// + Kx(r)+^ |^(r,t)p 



$(r,t), 



(5) 



which can be used to discuss the excitations from the equilibrium condensate wavefunction 
*o(r). 

To quantize a hydro dynamic-type description, we first replace the dynamical variables in 
(H) by the corresponding quantum mechanical operators. The density and velocity operators 



satisfy the exact commutation rule [|K| 



Ti 



(r),v(r')] =i— V5(r-r' 



m 



(6) 



This can be used in the analysis of the T = GP expression (f|) because almost all the 
atoms are in the condensate and hence p c (r) ~ p(r). (By the same token, it is difficult 
to quantize a hydro dynamic-type theory at T ^ when there is significant depletion.) In 
terms of the phase operator defined by (H), @ is equivalent to |T(J 

[p(r)J(v')]=i5(r-v'). (7) 

Using the operator version of the Hamiltonian in ([3]), it is straightforward to work out the 

Heisenberg equations of motion for both p and v making use of (H). One obtains [|9] 

dp(r,t) 



dt 



-V- [p(r,t)v(r,t)\ 



and 



m- 



<9v(r, t) 
dt 



-V 



V m (r) - p + g p{r,t) 



^vyp(r,t) i 

+ -mv (r, rj 



(8) 



(9) 



2m^Jp(r,t) ^ 

Assuming the validity of the Gross-Pitaevskii Hamiltonian given by (jl]), @ and @ are 
exact and can be used to discuss the dynamics of the condensate. They are equivalent to 
the GP equation in (|5|). 

It is useful to separate various operators into equilibrium and small fluctuation parts 

$(r,t) = $ (r) + 5$(r,t) 
p(r,t) = p (r) + 5p(r,t) 
0(r,t) = o (r)+50(r,t) 

v(r,t) =v (r) + <fv(r,*), (10) 



where the static condensate density profile is po( r ) = | < ^ ) o( r )| 2 and v (r) = hV4>o(r)/m. 
Before preceeding, we recall that the equilibrium values are given by the solutions of the 
time-independent GP equation [see (0)], 

h 2 v 2 



2m 



fM + V ex (r) + g p (r) 



$ (r) = 0, 



(11) 



where $o( r ) = |$o( r )l e ^ o( ' r ' ) - The static solution of (|8[) must satisfy 



V 



l*o(r)| 2 v (r) 



(12) 



which is equivalent to 



2V|$ (r)| ■ V0 o (r) + |$o(r)|V 2 <Mr) = 0. 



(13) 



Using this in ((TTJ) , the latter equation can be reduced to an equation for the amplitude 
l*o(r)| = 



h 2 V 2 1 

- /i + Kx(r) + ^oPo(r) + 2 mv o( r ) 



|$ (r)| = 0. 



(14) 



Solutions of (|i~3|) and ( |14|) with a spatially-varying phase 0o( r ) correspond to states with a 
non-zero superfluid velocity v (r) (such as vortices). 

In this paper, we are interested in small fluctuations of the static condensate $o( r )- I n 
terms of ([|a|) and (|10l), we obtain 



8p(r, t) = $o(r)5l>(r, t) + $o(r)^ f (r, t). 



(15) 



Writing $(r, t) = $o( r )[l + SA(r, t)]e t5 ^ r ^ in terms of Hermitian amplitude and phase 
fluctuation operators, we find 



5$(r,t) = $ (r)[5i(r,t) + iS}(r,t)] 



(16) 



and hence the density fluctuation operator in (|T3p is related (as expected) to the amplitude 
fluctuation operator 



Sp(r,t) = 2|$ (r)| 2 5i(r,t) = 2p (r)<L4(r, t). 



(17) 
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We have used the fact that 5 A = S A* and 5(f) = 5ft. With (0), the commutator in (0) can 
be rewritten in terms of fluctuation operators 

[5p(r),5$(r')]=i5(r-r f ). (18) 

The key simplification to be used in analyzing the Eqs. @ and (§) is to neglect the 
kinetic energy term associated with density fluctuations fi 2 V 2 \fp / "Imyfp in determining 
both the equilibrium and time-dependent solutions. This is referred to as the Thomas- 
Fermi approximation (TFA) [|13|,(J]. In the TFA, for v (r) = 0, the static density profile 
Po( r ) is related to the external trapping potential by fi£B 



[i = #oPo(r) + Kx(r), (19) 

where the chemical potential p is fixed by the total number of atoms N = J drpo(r) and 
the condition po( r ) > 0. This should be a good description if the number of atoms is large 
enough. More specifically, it is valid for QT3J Na/a HO 1, where ano = (h/muj )^ is the 
characteristic harmonic oscillator length for an external potential V eyL (r) = ^muj 2 r' 2 . One 
can show that this implies \i ^> huj . The use of the TFA precludes the consideration of a 
vortex state, where vq(t) = nh/mr± (r± is the distance from the axis of the vortex and n is 
an integer). In the presence of vortices, the kinetic energy terms in (fH| ) plays a crucial role 
in the determination of the spatial variation of |$o( r )| [0 • 

In the same approximation, (§) and (||) give (to first order in the perturbed quantities) 

dSp(r,t) 



dt 

and 



-V-[ Po (v)5v(v,t)] (20) 



m^^ = -W[g 5p(v,t)]. (21) 



To be explicit, we have neglected the following terms in the square bracket in 

V 2 ^o(ri 



AmJpo(r) [ Po(r) 



-5p(r,t) + V z 



Sp{r, t) 



Po{r) 



(22) 



relative to go5p(r,t). Combining (|20|) and fl21~|) gives 



dt 2 



and the equivalent 



<9 2 <50(r,t) 



%o(r)V5p(r,t) 
m 



%(r)V50(r,t) 

m 



(23) 



The solutions of ( f23| ) or (|24|) give the low-frequency condensate collective modes 
inhomogeneous Bose gas with a local condensate density Po( r ) given by (|T9|). 
In this approximate treatment, the GP Hamiltonian (Q) reduces to 



H = H + ^ I dr mp (r)5v 2 (r) + g^p 2 ^) 



where the ground state is described by 

H o = ~\j dr 9opl(r). 



(24) 



of an 



(25) 



(26) 



An expression identical to (^) is derived in § 5.2 of Ref. |J] by a different procedure. Clearly 
the quadratic form (^) can be diagonalized by expanding the operators Sp(r) and <50(r) [or 
<5v(r)] 



j 

*)^[B^(r)« i + B^(r 



(27) 



in terms of creation and annihilation operators aj and which create and destroy excita- 
tions of the condensate with energy hujj. These operators satisfy the usual Bose commutation 
relations: [d^d/] = [dj,dj,] = and [dj,dt,] = Sjji. The hermiticity of 5p and 5<p is en- 
sured by the form used in ([27]). The eigenfunctions ^-(r) are assumed to satisfy both the 
normalization relation 



dr ^*(r)^,(r) = 8 k 



3 * 



(28) 



and the completeness relation 



H^(r)^(r) =5(r-r'). 



(29) 



Using (p7[) with the help of (p9|), the commutator ([18]) can be shown to require that A*B 3 



-i/2. 



Assuming that ( |2"5| ) has been diagonalized to give 

H = H Q + Y. huJ J ( & j & J + I) > ( 3 °) 
the time-dependent operators in Heisenberg representation are then immediately given by 

dp(r,t) = [A^^e-^&j + A*V*(r)e^*a]" 
i 

50(r, t) = £ [s^We-^'dj + B;V$(r)e**'a}] . (31) 
j 

Using (|3~TD in immediately gives a second relation between the coefficients in (p7|), 
namely (for cjj > 0) 



ifiujjBj = goAj. 



(32) 



Using (|23|) or (|24D in connection with the expression (|3ll) , we see that the eigenfunctions 
ipj{r) and eigenvalues cjj are determined by 



V 



m 



i Po(r)V^-(r) 



Combining (|32|) with the relation Ajj-Bj = — i/2 leads to (for ujj > 0) 



Ihujj 

2^ 



' go 
2^,0;^ ' 



(33) 



(34) 



where, for convenience, A, is chosen to be purely imaginary and Bj is chosen to be real. 
Making use of the results (|34|) in fl2"7|), one can now verify after some algebra that the 



Hamiltonian fl25|) does indeed reduce to (|30|), as assumed. The harmonic Hamiltonian ( |30] ) 
describes the fluctuations of the condensate $o( r ) in terms of a non- interacting gas of Bose 
excitations. 

Eq. (^) has been solved by Stringari |J for both isotropic and anisotropic parabolic 
traps. For an isotropic harmonic potential, the solutions are (j represents the usual quantum 
numbers (n, £, m) for a spherical potential) 
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V»i(r) = CiP? n) (^) / F, m (0, 0) 6(i? - r), (35) 
where 0(sc) is the step function. The associated energy eigenvalues are found to be 

Uj = uoq (2n 2 + 2ni + 3n + i) 1 . (36) 

n 

In ([55]) , Yg m (9,(f)) are spherical harmonics, while P/ 2 ™^) = «2m(^ 2m is a polynomial 

m=0 

with coefficients satisfying the recurrence relation: a2 m +2 — —&2m{ n ~ m)(2£ + 2m + 2n + 
3)/(m + 1)(2£ + 2m + 3), with ao = 1. In the Thomas-Fermi approximation (|T9|), the 
chemical potential is given by p = ^mu 2 R 2 , where R is the radius at which po( r ) vanishes. 



This implies that the integration in fl28|) is over a sphere of radius R. One also finds directly 
that p = 15ah 2 N j '2mR? , which can be used to give R in terms of N and a. One can show 
that the TF hydrodynamic approximation should be valid for u 3 < p || . The normalization 
factor Cj is determined via ( p8|) , 



1 _ . _ r ,„ ^ -, 2 ~| 2 



c, = c n , m = \ R 2 ^ A I dx x 2 ^ 2 Pt\x) ~\ . (37) 



For n = 0, one has c olm = [{21 + 3)/iF +3 ]i 

For an isotropic parabolic trap, fl35|) gives the normal-mode basis functions which appear 
in 5p(r,t) and <50(r,t) in (^). The normal- mode frequencies (^) are independent of the 
interaction strength because in ([£]), we have used the relation g Q p (r) = p — V cx (r) given by 
(P^|). The corresponding results of (^) and flSB] ) for a uniform Bose fluid |TD[ are ^-(r) ~~ > 
^k(r) = ^e* k ' r and — > u(k) = ck with c 2 = g p /m. In this case, (|19|) gives the usual 
Bogoliubov result p = goPo- We note that only finite-energy solutions of ( |3~3"|) are involved 
in the fluctuations ||. A zero energy solution would correspond to the ground state value. 
The n = 0, £ = solution of ([35]), which corresponds to V'jl 1 ') = constant and LUj = 0, is 
excluded in all summations such as (|27|) . For further discussion of the zero energy solution, 
see the end of Section V. 
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III. GREEN'S FUNCTIONS AND DEPLETION OF CONDENSATE 



Writing the quantum field operators (see Ch. 3 of Ref. [O) in terms of condensate 
and non-condensate contributions ^(r) = $o( r ) + 4>( r ), the diagonal single-particle Green's 
function Gn is given by [using (|i5|) -([i7|)1 



= -(5<f> t (r,i)5<!>(r , )> 



$S(r)$o(r) 



{SA\r,t)SA(r')) + (<ty t (r,t)<ty(r')> 



+ i(5A^(v,t)5(j)(r')) - i(5ft(r,t)5A(r')) 



(38) 



Using the normal-mode expansion fl27|) , we find after some algebra 



Gn(ry,t) = -$S(r)* (r')E 
1 



+ 



go 



+ 



Po 



Po(r)po(r') 2^- 
(r) p (r')J J 



(39) 



where we have defined 



C j± (v,r',t) =V;(r)^(0^(^)c^'±Vi(r)V?('') ^Vj) + 1 



(40) 



and N°(ujj) = [exp(f3hu>j) — l]^ 1 is the Bose distribution function, with /3 = 1/ksT. In (|39"1), 
the first term involves amplitude fluctuations, the second term involves phase fluctuations, 
while the third and fourth terms involve the coupling between the amplitude and phase 
fluctuations. Strictly speaking, this result is only valid at T ~ 0, since we started with the 
Gross-Pitaevskii description (ffl) which assumes all the atoms are in the Bose-condensate. 
The time-ordered Green's function equivalent of ( j39"D is easily worked out. 
In an analogous way, one can calculate G22 to be 



= -<&o(r)$*(r')£ { 



+ 



0o 



1 1 

+ 



Mr) M r 



80o Mr)po(r') 2huJ j 
C,_(r,r',t)j, 



(41) 
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and the off-diagonal single-particle Beliaev Green's functions are 



Gi 2 (r,r',t) 



-<<5$ t (r,t)5<i> t (r , )> 
-*o(r)<&o(r')£ 



hujj 



p (r) po(r') 



%o Po(r)p (r') 



9o 
2hujj 



C i+ (r,r',t) 



(42) 



and 



Gf 2 i(r,r , ,t) = -( < J$(r,t) < yd(r')) 
= -$ (r)$ (r')E^ 



9o 



+ 



. 8^0 Po(r)p (r') 2hujj_ 
Q-(r,r',t)} 



C,- + (r,r',t) 



(43) 



.Mr J Pol 

As an application of the above results, we calculate the depletion of the condensate 
density near T = 0. Using (|39|) and ([40]), one obtains 



(p(r)) = -G n {r,r,t = 0) 



5>; 



:r)i 2 i 



./ , 9oPo(r) 



(T ~ 0). 



(44) 



_4g p (r) hujj 

For a spatially uniform Bose-condensed gas (po(r) = constant and c^- — > u>(k) = ck), the 
phase fluctuations (second term) are the dominant contribution to the expression in (|44"|). 
It then reduces to [fLOU 



1 f dk g p 
P = P ~ Po = - 



(45) 



2 J (2vr) 3 ^(k) 

In contrast, the phase fluctuations are not so dominant in a trapped Bose gas because their 
energy is finite. We note that (34) is similar to recent results based on a variational solution 



of the full Bogoliubov equations [16 



As discussed by Stringari ||, the collective mode energies in ( plop in the n = case 
reduce to the simple expression Uj — > ojoe = uo\^i, for i finite. The ^-channel contribution 
in (54) for these n = "surface modes" || is given by 
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'VI a 2 „ n . 1 - / 



= + 1)(2£ + 3 > •* [t^W + w£ " 2 1 • (46) 



where 2 m \Yim(0, 0)| 2 = (2£ + l)/47r has been used and we have defined r = r/R, auo = 
ano/R- The n = 1, £ = contribution to the depletion is (<ui,q = v^5^o) 



647T.R 3 V 3 / V 1 - f 2 Vla 2 HO 



For illustration, in Fig. |1] we show the depletion (p(r)) as a function of r from the lowest- 
energy modes described by ( fjE|) and ( f4"T| ) [ |H| ]. We have assumed iV = 10 5 trapped atoms, 
an s-wave scattering length a = 50A and auo — 10 4 A. These parameters give R = 6x 10 4 A 
14]. The local condensate density is given by (19), namely 

PoW = — ^-(l- r 2 ) 6(1- f). (48) 

For the parameters used in Fig. [I], we find po(r) = (1 — ^ 2 ) 2.8 x 10 14 cm -3 and thus the 
T = condensate depletion is very small, as expected. For these parameters, we note that 
fi = | (^~~) h^o — 18/itUo an d hence these frequencies satisfy the condition u> < \x for which 
the TF approximation is valid l|. 



IV. DYNAMIC STRUCTURE FACTOR 

We next calculate the light scattering cross section which is proportional to the dynamic 
structure factor S(q,u). The latter is given by the Fourier transform of the density-density 
correlation function. We recall that the density can be split into contribution involving 
atoms being excited in/out of the condensate and contribution from non-condensate atoms. 
Thus we have (see Ch. 3 of |TT| ) 

p(r) = $(r)ip(r) 

= |$ (r)| 2 + [^(r)^(r) + $o(r)^(r)] + ft(r)$(r) 

= p (r) + 5p c (r)+p(r), (49) 
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where p(r) = ij)*(r)ij)(r) is the non-condensate local density operator. In the present dis- 
cussion of a dilute Bose gas, it will be assumed that one can neglect p(r) as being small at 
T = 0. The density- density correlation function due to the condensate density fluctuations 
is easily calculated to be [using (|17D, (p8[), and (|3^)] 



(5p c (v,t)5p c (r')) = 4p (r)po(r / )(5i(r,t)5i(r')) 
= £|^+(r,r',t), 

j 290 



(50) 



where the function Cj + (r, r',t) has been defined in (0). Thus for the problem being con- 
sidered, the condensate contribution to >S(q, uj) is given by 



S c (q, uo) (x J dre 

,TlUJA 



dte wt {6p c (T,t)5p c (T')) 



(51) 



where 



dre- iq >,-(r). 



(52) 



This is for a momentum transfer ftq and energy transfer /icu given to the Bose condensate. 

Since we are only dealing with condensate fluctuations (density fluctuations involving 
atoms entering/leaving the condensate), it is trivial that the poles of 5 c (q, uj) in (^T|) are 
identical to those of the single-particle Green's functions in (|38|) and (0). However, this 
sharing of poles is a general feature of Bose-condensed fluids even when one includes non- 



condensate contributions (see Ch. 5 of Ref. 11 ). While (pTD is based on a calculation which 



is only valid at T ~ 0, we expect that more generally, it will approximately describe the 
contribution related to exciting excitations out of the condensate. In this context, it is useful 
to give the two limiting expressions 



S c (q,uj) 



OJ i 



if ksT <C hujj 



£^%i(q)l 2 %- 
i 290 

— lV ; i(q)| 2 [^( w - Wj) + 5(uu + ujj)} if k B T > hujj 
i 2 9o 



(53) 
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Of course, ([51]) and (|53|) do not include the scattering from thermally excited atoms present 



at T ^ or the two-excitation ( "multiphonon" ) contributions flTT| , |T2|| . These are briefly 
discussed at the end of this section. 

In the analogous calculation of 5* c (q, to) for a weakly interacting homogeneous Bose- 
condensed gas OJT2I], one picks up a single sharp peak at the quasiparticle energy uj(q) cor- 



responding to the momentum transfer involved. In a trapped Bose gas, in contrast, S^q, uj) 
in is seen to be a weighted sum of all the normal modes describing the fluctuating 
inhomogeneous condensate, of frequency ujj and with a weight proportional to ujj\ipj(q)\ 2 . 
Using the well-known expansion (for q parallel to z-axis) 



e - H , = ^ ^ 47r(2£ + l){ _iy h{qr)Ym (e), (54) 



where ji is the f-order Bessel function, one can use (|35|) to directly obtain the Fourier 
transform of the eigenfunction ipj{r), namely 

VWm(q) = c nem (-iY5 m>0S /47r(2£+l) R e+3 f dx x l+2 P^ n \x)j,{qRx), (55) 

•> 

where the normalization factor c n £ m is given by (|37|). We note that as a result of (|54]), only 
the states V'i( r ) with m = have a finite weight in S" c (q, uj) in (|5"T|). 

In Fig. we have plotted uij\ijjj(q)\ 2 which appears in the formula (|5l|) for 5* c (q, uj) as a 
function of the dimensionless wavevector qR for excitations (n, £) = (0, 1), (0, 2), (0, 3) and 
(1,0), with m = in all cases. As one can see from Fig. 0, the strongest weights for these 
collective modes appear for qR < 6. As an example, for the parameters used in Fig. [1] to 
give R = 6 x 10 _4 cm, this means that ideally, the momentum transfer q in a light-scattering 
experiment should be q ^ 10 4 cm _1 in order to pick up the strong spectral weight from 
the low-energy collective modes. For qR > 1, the density fluctuations being probed have a 
wavelength smaller than the size of the condensate 2R. The Thomas-Fermi approximation 



(|j~9|) we have used should be adequate. In contrast, the TF approximation is probably not 
very good for qR -C 1 since the results could be sensitive to the abrupt vanishing of po( r ) 
for r > R. 
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In Fig. ^ we plot the T = dynamic structure factor S" c (q, u) given by (|53|). For clarity 
and as a description of finite energy resolution, we have broadened the delta function S(uj) 
using a Lorentzian, with a width of T = 0.05o;o- The momentum transfer corresponds to 
qR = 2. In Fig. |], we give results for a higher momentum transfer qR = 10 using the same 
arbitrary units as in Fig. |3], showing a much weaker light-scattering intensity. 

Besides the condensate contribution to S"(q, uj) given by (|51~D , we can use our hydrody- 
namic approximation to estimate the contribution of the non-condensate part. Using (|49|). 
the density-density correlation function arising from the non-condensate fluctuations can be 
expressed in terms of the single-particle Green's function G a p discussed in Section III, 

(p(r,t)p(r')> = (5&{r,t)5${r,t)5&(r')5$(r')) 

= (5$ t (r,t)5$(r / ))(5$(r,0^ t (r / )> + (5$(r,t)5$(r'))(5& (r,t)5& (r')) 

= G n (r, r\ t)G 22 (r, r', t) + G 21 (r, r', t)G 12 (r, r', t). (56) 

After some algebra, using (|3^)-(|43|), we obtain the T = contribution due to creating two 
excitations, 

S(q,to) oc / dre- iq - r / dr'e^' f°° dte tujt (p(r, t)p(r')) 

J J J — OO 

= 5321^(0)1^(0;-^-^), (57) 



where we have defined 



%(q) = / rfre " 



iq-r 



ibJrUJr). (58) 
This result is the generalization of a similar calculation for a uniform Bose gas given in 



Ref. [11]. Of course, at T = 0, this contribution is much smaller than the condensate 



contribution S c (q,u) in (^T|) since most of the atoms are in the condensate, i.e., p p c . 



Using (p6|), one can see that S(q, o;) contributes only when u > 2u . 

The general expressions ( |5T| ) and ( |5"7D given in this section are also valid for anisotropic 
parabolic wells. Stringari has discussed the solutions ipj(r) of ( p3|) for traps with axial 
symmetry about the z-axis. In this case, the azimuthal quantum number m still characterizes 
the normal-mode eigenvectors if}j(r). 
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V. RELATION TO BOGOLIUBOV COUPLED EQUATIONS 



The standard approach to study the dynamics of a T = Bose condensate (uniform 
or non-uniform) is based on the Bogoliubov approximation 0,B]. The non-condensate field 
operators are decomposed into positive- and negative-energy Bose fluctuations 



O i<i>o(r) 



EE 



■iEjt/h 



a , 



(59) 



where aj, d] satisfy Bose commutation relations. The amplitude functions Uj(r) and Vj(r) 
in (|59|) are given by the solutions of the coupled eigenvalue equations 



L Uj (r) - g \%(r)\ 2 Vj (r) = E jUj (r) 
L* Vj (v) - g \%(r)\ 2 Uj (r) = -E jVj {v 



(60) 



with 



L = - — [-ihV + mv (r)f -fjt + V ex (r) + 2g \%(r)\ 2 . 
2m 



(61) 



As noted by Fetter |I6[, the solutions of (|60| ) for Uj and Vj are basically linear combinations 
of 5p and 8cj) in the TF approximation studied in this paper. Within this approximation and 
setting v (r) = 0, ( JIB] ) and ([H]) give 

~6p(r,t) 



i)(r,t) =5§(r,t) = $ (r 
Comparing terms in (1591) with (pTj), (I 



2po(r) 
immediately gives 



+ i6(/>(r,t) 



(62) 



^•(r) 



ffoPo( r ) 
\\ 2huj j 



+ 



hujj 



%oPo(r) 



ffoPoW 
\ 2^,0;^ 



8fi'oPo(r) 



^i(r), 



(63) 



with Ej = hujj. Here we have used |$o( r )| = yPo( r )- One can easily verify using ( |28"D and 
(fZ9|) that Uj(r) and fj(r) given by ( pB| ) satisfy the usual orthogonality and completeness 
relations. 
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For completeness, it is useful to show explicitly how the Bogoliubov equations (|60| ) can 
be used to derive the key equations fl20| ) and fl2~T|) on which the TF hydrodynamic description 
is based. Using ( |T9| ) , the operator in (|6T| ) can be written as L = Lqp + go\&o( r )\ 2 , where 

h 2 v 2 



L GP = 



2m 



p + V ex (r)+g \<f> (r)\ 2 . 



(64) 



Subtracting and adding the two equations in (^) gives 

Lgp [uj(r) - Vj(r)] + 2g oPo (r) [uj(r) - Vj (r)] = Ej [ Uj (r) + Vj (r)} 
L GP [ Uj {r) + Vj(r)} = Ej [uj(r) - ^(r)] . 



(65) 
(66) 



Making use of (|63|), d27|) , and fl34|) , one finds by direct comparison 

5pj{v) 



l*o(r)| 



Uj(r) + Vj(r) = 2i, 



floM r J 



^(r)=2i|$ (r)|^(p) 



(67) 



where, using (|27j) and fl34|), we have defined 5pj(r) = ^-^(r) and <5</>j(r) = ^^-^(r). 
Taking into account fll9|), we see that the first term on the r.h.s. of (^) reduces to 



2m 



V 2 [5p i (r)/|$ (r)|]. Since this is omitted in the TFA [see (f22|)], gives 

iEjSfair) = g 5pj(r). 



(68) 



Using this in (|31~D , this is equivalent to (|2l| ) or 



<9i 

Similarly, substituting ( |B~7D into (|66|) gives 



-9oSp(r,t). 



(69) 



$ (r)|V^[|$o(r)|^ i (r)]=i^p i (r). 



(70) 



Within the TF approximation (i.e., neglecting terms involving V 2 |$o( r )l)) d^0|) is equivalent 
to 



\M*)\ 2 to i {T)]=i^5p j {T) 



(71) 
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which, using (pi]), is equivalent to (pC|). 

We note that while and fl66|) describe fluctuations of the equilibrium condensate, 
formally they also exhibit a special solution (labelled by s) corresponding to u s (r) = v s (r) 
with zero energy E s = If one identifies u s (r) with |$(r)|, fl66|) is then equivalent 

to the static GP equation in This zero-energy solution still appears in the TF hy- 

drodynamic approximation described by (|67|), in which case we see that Sp s (r) = and 
u s (r) = i|$ o ( r )|^0s( r )- This is a zero energy solution of (|24l) corresponding to a phase 
change 50 s (r) which is independent of r. This zero frequency fluctuation of the order pa- 
rameter has a very simple physical meaning, namely it involves uniform (rigid) phase change 
5(f> s at all points, but with no associated amplitude (or density) fluctuation (5p s (r) = 0). 
This is the expected zero frequency mode in a system with a Bose broken symmetry. 

VI. CONCLUDING REMARKS 

The quantized expressions of density and phase fluctuations give a simple way to calculate 
the contribution of low-energy condensate fluctuations to various observable quantities. The 
results given in this paper are based on the "hydrodynamic description" in terms of the 
amplitude and phase of the condensate wavefunction at zero temperature, as formulated 
by Stringari || for trapped Bose fluids, in which the kinetic energy contribution due to 
density fluctuations is neglected. While this Thomas-Fermi hydrodynamic description is 
only adequate for the low-frequency modes (ujj < /x), it is much simpler than dealing with 
the full set of coupled Bogoliubov equations given by (|60|) . 

Calculations based on these normal modes allow one to exhibit the role of amplitude and 
phase fluctuations of the condensate in a very explicit fashion, complementing a more com- 
plete calculation based on numerical solutions of the Bogoliubov equations [§-[§]. We have 
illustrated this formalism by evaluating the single-particle diagonal and off-diagonal Green's 
functions, the T = condensate depletion, and the dynamic structure factor S'(q, o>), taking 
into account the contribution of the lowest-frequency modes. These specific calculations are 
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carried out for an isotropic parabolic trap but can be generalized for an anisotropic trap 
using the results in Ref. 0. 
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FIGURES 




FIG. 1. Plot of the depletion given by (^J) of the local condensate density due to the con- 



tributions of several low-energy excitations specified by (n,£) 17]. The parameters used are: 
a HO = 10 4 A, a = 50A, and N = 10 5 , which give R = 6 x 10 4 A. 
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FIG. 2. Plot of Uj | ipj (q) | 2 given in (^) as a function of dimensionless wavevector qR, where 
R is the size of the condensate. These functions determine the weight of the light-scattering cross 
section in ( |53| ) of the corresponding collective modes of energy uj n £o (see Figs. || and Q). 
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FIG. 3. Plot of the T = dynamic structure factor 5(q, u>) vs to, using a Lorentzian broadening 
of the delta function peaks (T = 0.05a>o) and a momentum transfer qR = 2. Only the contributions 
from the low-frequency modes (n,£, 0) are included. 
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FIG. 4. Same plot as Fig. ^, for qR = 10. Note the much smaller intensity compared to the 
results for qR = 2 in Fig. ||. 
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